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Abstract 

Ordinary-sterile nentrino oscillations can generate significant nentrino asymmetry in 
the early Universe. In this paper we extend this work by compnting the evolution of 
neutrino asymmetries and light element abundances during the Big Bang Nucleosyn¬ 
thesis (BBN) epoch. We show that a significant electron-neutrino asymmetry can be 
generated in a way that is approximately independent of the oscillation parameters hirP 
and sin^ 29 for a range of parameters in an interesting class of models. The numerical 
value of the asymmetry leads to the prediction that the effective number of neutrino 
flavours during BBN is either about 2.5 or 3.4 depending on the sign of the asymmetry. 
Interestingly, one class of primordial deuterium abundance data favours an effective 
number of neutrino flavours during the epoch of BBN of less than 3. 


I. Introduction 


The possible existence of sterile neutrinos can be motivated by the solar neutrino, 
atmospheric neutrino and LSND experiments]^]. There are also interesting theoretical 
motivations for the existence of light sterile neutrinos. For example, if nature respects 
an exact unbroken parity symmetry, then three necessarily light mirror neutrinos must 
exist 0. In view of this, it is interesting to study the implications of ordinary - sterile 
neutrino oscillations for both particle physics and cosmology. The effects of ordinary - 
sterile neutrino oscillations in the early Universe are actually quite remarkable. It turns 
out that for a wide range of parameters, ordinary - sterile neutrino oscillations generate 
a large neutrino asymmetry [Q (see also |^). (A large neutrino asymmetry implies that 
the universe has a net nonzero lepton number given that the electron asymmetry is 
necessarily small due to the charge neutrality requirement.) One important implication 
of this result is that the bounds on ordinary - sterile oscillation parameters that can 
be derived mainly from energy density considerations during Big Bang Nucleosynthesis 
(BBN) are severely affected (see Ref.for a detailed analysis). However, electron 
lepton number can also affect BBN directly through the modification of nuclear reaction 
rates. It is this issue that we will study in this paper. 

In a previous paper |^, we showed that for a wide range of parameters, the evo¬ 
lution of lepton number can be approximately described by a relatively simple first 
order integro-differential equation. We called the approximation used there the “static 
approximation” because it holds provided that the system is sufficiently smooth. The 
static approximation is valid in the region where the evolution of lepton number is 
dominated by collisions. In particular, for the temperature at which lepton number is 
initially produced, this approximation is generally valid for \5m?\ ~ 10“^ eU^0. How¬ 
ever, it is not expected to be valid for temperatures much less than the temperature at 
which lepton number is initially generated. This is because the static approximation 
discussed in Ref.[Q does not incorporate the MSW effect 0, which is in fact the dom¬ 
inant process affecting the evolution of lepton number at low temperatures. For the 
application considered in Ref.^, the evolution of lepton number at low temperatures 
was not required. However, for the application of the present paper the accurate evo¬ 
lution of lepton number to temperatures T ~ 1 MeV is necessary in order to study its 
precise effect on BBN reaction rates. 

The outline of this paper is as follows. In section H we set the scene with a brief 
review of the effects of neutrino asymmetry on BBN. In section HI we develop a simple 
formalism describing the evolution of lepton number at low temperatures where the 
MSW effect is important. This work can also be viewed as an extension of our previous 
study P, where the evolution of lepton number at higher temperature was studied in 
detail. In this section we also examine the implications for BBN of direct electron asym¬ 
metry generation by Uf. — Ug oscillations. In section IV we examine a more interesting 
scenario where electron asymmetry is generated indirectly. In section V, we provide a 
check of our simple formalism (of section HI) by numerically solving the exact quantum 
kinetic equations. Finally in section VI we conclude. 
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II Electron neutrino asymmetry and BBN 


Standard BBN can give a prediction for the effective nnmber of nentrino flavonrs, 
, present dnring nncleosynthesis. This prediction depends on the baryon to photon 
nnmber-density ratio, rj, and the primordial helinm mass fraction, Yp. A precise deter¬ 
mination of the primordial denterinm abnndance will provide a qnite sensitive measnre- 
ment of rj. Once rj is known, the effective nnmber of nentrinos present dnring nncleosyn¬ 
thesis depends only on Yp. At present there are two conflicting denterinm observations 
in different high-redshift low-metallicity qnasistellar object absorbers. There is the high 
denterinm resnlt of Ref.|3 which snggests that D/H = (1.9±0.4) x 10“"^. On the other 
hand, in Ref.[|| the low denterinm resnlt of D/H = [2.3 ± 0.3{stat) ± 0.3{syst)] x 10“^ 
is obtained. The implications of these results for the prediction of have been 

discussed in a number of recent papersj^. Depending on which of these two values 
of the deuterium abundance is assumed, different predictions for y are obtained. The 
high deuterium result leads to r; ~ 2 x 10“^°, while the low deuterium result leads to 
r] ~ 7 X 10“^° [0. Each of these predictions for rj, together with the inferred primordial 
abundance of '^He, allows a prediction for to be made0. According to Ref.||TD| 
for example, the high deuterium case leads to 


= 2.9 ±0.3, 


( 1 ) 


while the low deuterium case leads to 


= 1.9 ±0.3, 


( 2 ) 


where the errors are at 68% C.L. The minimal standard model of course predicts that 
= 3. Thus, if the low deuterium result were correct, then new physics would 
presumably be required]^. Of course, estimating the primordial element abundances 
is difficult and it is possible that the primordial helium abundance has been underesti¬ 
mated (in otherwords, even if the low deuterium measurement is correct = 3 is not 
inconsistent). Fortunately the situation is continually improving as more observations 
and analyses are done. In the interim it is useful to identify and study the types of 
particle physics that can lead to < 3. 

One possibility is that the electron lepton number is large enough to signihcantly 
affect BBN (i.e. ~ 0 .01)1^. The relationship between an electron neutrino asym¬ 

metry and the effective number of neutrino species arises as follows. A nonzero electron 
neutrino asymmetry modihes the nucleon reaction rates {n+Vf. ^ p±e“, n+e^ ^ P+^e) 
which keep the neutrons and protons in thermal equilibrium down to temperatures of 
about 0.7 MeV. A modihcation of these rates affects the ratio of neutrons/protons 
and hence changes the prediction of Yp[|^]. A change of Yp can be equivalently ex¬ 
pressed as a change in the effective number of neutrino species, , present during 

nucleosynthesis. These quantities are related by the equation (see e.g. Ref.|1IB|) 


5Yp ~ 0.0125Ar5^. 


( 3 ) 
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The effect of the electron neutrino asymmetry on the primordial helium abundance is 
most important in the temperature region 

0.4 MeV ~ T ~ 1.5 MeV (4) 

where the reactions n + i/g ^ P + and n + p + z/g fix the neutron/proton ratio. 
For temperatures less than about 0.4 MeV, these reaction rates become so slow that 
the dominant process affecting the neutron/proton ratio is neutron decay. Note that 
the Helium mass fraction Yp satisfies the differential equation ||14||, 

-^ =-X(n ^ p)Yp +X{p ^ n){2-Yp), (5) 

where A(n —> p) [A(p —> n,)] is the rate at which neutrons are converted into protons 
[protons are converted to neutrons]. For temperatures in the range of Eq.(^), A(n —>■ 
p) ~ A(n + z/e —p + e~) + \{n + e"*" —>■ p + z/g) and A(p —n) ~ A(p + Ue ^ 
n + e~^) + A(p + e~ —> n + z/g). The reaction rates (per nucleon) are obtained by 
integrating the square of the matrix element weighted by the available phase-space. 
For example, the rate for the process n + z/g —>■ p -1- e“ is given by 

\{n+v^ p+e-) = j 

( 6 ) 

where fi{Ei) is the Fermi-Dirac distribution fi{Ei) = [exp{Ei/T) -|-1] These reaction 
rates are modified in the presence of significant electron neutrino asymmetry. If the 
neutrino asymmetry is produced at temperatures above about 1.5 MeV and is constant 
over the temperature range of Eq.(^ then we only need to add in the appropriate 
chemical potentials pu and pp to the distributions fip and fp. 

Ill Neutrino oscillation generated neutrino asymmetry 

We now discuss the effects of neutrino oscillations, assuming that a sterile neutrino 
exists. Our convention for the neutrino oscillation parameters, and sin^ 29o, is as 
follows. For Ua — oscillations (with a = e,p,T), the weak eigenstates z/^ and z/^ are 
linear combinations of mass eigenstates Va and 

Ua = cos OoUa + siu OqUi,, Z/^ = — siu + COS 6*oZ/fe, (7) 

where the vacuum mixing angle 6q is defined so that cos 26^0 > 0. Further, we define 
the oscillation parameter by = ml — ml- Also, the term “neutrino” will 

sometimes include anti-neutrino as well. We hope that the correct meaning will be clear 
from context. 

Ordinary - sterile neutrino oscillations can create significant lepton number provided 
that < 0 and ~ eV^. For full details, see Refs.[^, ||, §|. In the following 

we consider z/q — Ug oscillations in isolation. It is important to note that this is not in 
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general valid because the effective potential (see below) depends on all of the lepton 
asymmetries. However, it is approximately valid for the ordinary - sterile neutrino 
oscillations which have the largest 

The effective potential describing the coherent forward scattering of neutrinos of 
momentum p = |p| ~ T' with the background is||I^, ^ 


v; = V„{T,p, L‘»>) = + b), 

where the dimensionless functions a and b are given by 




( 8 ) 


(9) 


and (C(3) ~ 1.202 is the Riemann zeta function of 3, Gp is the Fermi constant, M-\v is 
the W-boson mass, A^, ~ 17 and ~ 4.9 [0. The quantity LA) is given by 


+ Li,^ + Ly + + p. 


( 10 ) 


where — np^)/n^ with n* being the number density of species i. In kinetic 

equilibrium n* and hence L^') is in general a function of the independent variables /ij (the 
chemical potential) and T. For the situation we will be considering, the asymmetry L^') 
quickly becomes independent of its initial value (see Ref.[Q for a complete discussion). 
This effectively means that /Xj is not an independent variable but rather it becomes a 
function of T. The asymmetry LA) is thus essentially a function of T only, and a and b 
are functions of p and T only. The quantity p~LAr/2isa small term (~ 10“^°) which 
arises from the asymmetries of baryons and electrons. The matter mixing angles are 
expressed in terms of the quantities a and b through ^ 


sin^ 29m. 
sin^ 26m 


sin2 20^(T,p,L(“)) 
sin2 20™(T,p,L(“)) 


sin^ 26*0 

sin^ 26*0 + {b — a — cos 26*o)^ ’ 
sin^ 26*0 

sin^ 26*0 + {b + a — cos 26*o)^ 


( 11 ) 


Note that the MSW resonance occurs for neutrinos of momentum p when 6m = 7 r /4 
and for antineutrinos of momentum p when dm = 7 r/ 4 , which from Eq. (|Tl|) implies that 
b — a = cos 26*0 and b + a = cos 26*o, respectively. 

If sin^ 200 -C 1, then it can be shown that oscillations with 6 < 1 create lepton 
number while the oscillations with 6 > 1 destroy lepton number P, |^. Since ( 6 ) ~ T®, 
it follows that at some point the lepton number creating oscillations dominate over 
the lepton number destroying oscillations (where the brackets (■ ■ ■) denote the thermal 
momentum average). We will call this temperature Tc. It is given roughly by the 
temperature where {b) = cos 20 o — 1 , that is 


Te ~ 13(16) 



MeV, 


( 12 ) 
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for Vf. — Ug oscillations. 

It is important to note that there are two distinct contribntions to the rate of 
change of lepton nnmber. One contribntion is from the oscillations between collisions. 
The other is from the collisions themselves which deplete nentrinos and anti-nentrinos 
at different rates in a CP asymmetric backgronnd. It tnrns ont that for T ~ Tc, 
lepton nnmber evolntion is dominated by collisions for the small vacnnm mixing angle 


case we are considering assnming that ~ 10 ^ eV'^ i I, S 0- Oscillations 

between collisions, and in particnlar MSW transitions, can be ignored for T ~ Tc 
becanse the interactions are so rapid that the mean distance between collisions L^nt 
is mnch smaller than the matter-oscillation length TJ^c (and conseqnently the nentrino 
cannot evolve coherently throngh the resonance). To see this note that the amplitnde of 
the oscillations at the MSW resonance is given ronghly by sin^ Lint/2L^g^, where Lint ~ 
1 /{GjpPresT'^) is the interaction length at the resonance and T™c ~ ‘^Pres /(sin 26*o|5m^^|) 
is the oscillation length at the resonance H. So, at the resonance. 


sin 200 15m; 




2 I 

as I 


G^PresT^ ‘2pr 


~ 10^ sin 200 


n 6 


(13) 


where we have nsed the approximation — (p) — 3.15T for the resonance momentnm. 
Thns sin^ Lint/‘2L'^g^ <C 1 for T ~ Tc, provided that sin^ 20o ~ 10“"^, and so the MSW 
transitions are heavily snppressed. However, ~ 1/T® rapidly increases as T 

becomes lower. Also, for T ~ Tc it tnrns ont that Pres/T ~ 0.8 (see later). Taking these 
factors into acconnt, Lint/L'^^^ ~ 1 for T ~ Tc/2, provided that sin^ 20o ~ 3 x 10“^°. 
In this case the MSW effect will not be snppressed if the oscillations are adiabatic. 
Fnrthermore the oscillations are generally adiabatic for the parameter space of interest 
in this paper (which tnrns ont to be |5m^^| ~ 10“^ eV"^). 

The key task in this paper is to analyse the evolntion of lepton nnmber for T ~ Tc/2 
when MSW transitions become important. (The effect of MSW transitions was noted in 
Ref, bnt since the evolntion of lepton nnmber to low temperatnres was not reqnired 
for the application considered there, we did not stndy the effect of MSW transitions in 
any detail, except to note that they keep growing like 1/T^ for low temperatnres). 
Now, throngh the qnantnm kinetic eqnations [^, the evolntion of lepton nnmber can 
be calcnlated in a close to exact manner, inclnding the effects of both collisions and 
oscillations between collisions. However, these complicated conpled eqnations have two 
notable drawbacks. First, they do not fnrnish as mnch physical insight as one might 
wish. Second, they are impractically complicated when one wishes to consider a system 
of more than two nentrino flavours. Since the physics of MSW transitions is the essence 
of how lepton number evolves during the BBN epoch, and since we will later need to 
consider a system of four neutrino flavours, we now pursue a very useful approximate 
approach instead of employing the full quantum kinetic equations. We will along the 
way check the veracity of our approximate approach by comparing results with those 
obtained from the quantum kinetic equations in the two flavour case (see section V). 
This will give us conhdence in the use of our appproximate formalism in the four-flavour 
case considered later in this paper. 
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For definiteness we will assume that the lepton number created at the temperature 
T = Tc is positive in sign In this case note that a,b > 0 given also that < 0. 
The momentum of the anti-neutrino oscillation resonance (obtained from the condition 
b + a = cos 26*0 — 1) typically moves to quite low values, Pres/T ~ 0.8 (for T ~ T^). 
In contrast, the neutrino oscillation resonance momentum obtained from —a -|- 6 ~ 1 
moves to a very high value, Pres/T ^ 1 (see Figure 2 in section V for an illustration 
of this). As 6 ~ (6) ~ T® becomes smaller, the neutrino momentum resonance Pres/T 
very quickly becomes so high that its effects can be neglected, because the resonance 
occurs in the tail of the neutrino momentum distribution. Thus, for T < T^, we can to 
a a good approximation ignore the neutrinos, and simply study the effects of the MSW 
transitions on the anti-neutrinos. In this case all the anti-neutrinos which pass through 
the resonance are converted into sterile neutrinos and vice-versa (the MSW effect). The 
rate of change of lepton number is thus related to the number of anti-neutrinos minus 
the number of sterile anti-neutrinos which pass through the resonance. Note that this 
rate is independent of the precise value of sin^ 26*o provided that sin^ 26*o <C 1. Under 
this assumption 


dLi. 


Np 


Np 


dT 


T 




d 

It 


Pr 


(14) 


where W describes the momentum distribution of species i, so that n* = Nidp. In 
thermal equilibrium. 


Np = 

Uot 


1 


p 


27r^ 1 - 1 - exp 


(15) 


In Eq. ([T^) the factor Td{pres/T)/dT ~ dpres/dT—p/T = dpj.es/dT—dp/dT is the rate at 
which Pres changes relative to the neutrino momentum (for neutrinos with momentum 
P ~ Pres)- Note that ~ 0 if the number-density of sterile neutrinos is negligible. 
The functions W in Eq.(|^ are evaluated at the resonance momentum, Pres-, obtained 
as a function of Ly^{T) and T from the resonance condition a ~ cos20o — 1, 

Pres[T, Ly^{T)] ^ g 

T 8Ci3)^/2GFTTj 

where we have considered the case p, Ly^ -C Ly^ for (3 ^ a. Note that this expression 
is only valid for T ~ Tc/2 where the 6-term can be neglected. Using 


Pr 


T 


d{pres/T) d{pres/T) d{pres/T) dLy^ 


dT 


dT 


+ 


dL„ 


dT 


= -4 


Pres Pres 


dL,. 


T2 TL. dT 


Eq.(|T4D yields 


dT], 


-AX Pres/T 


(17) 


(18) 


dT T + Xpres/LyJ 

where we have assumed d{pres/T) /dT < 0. The useful dimensionless quantity X is 
given by 


X ^ X[T,p,PpST),NpXT)] = - {Np^ - Np,) 


n.. 


(19) 
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and it is evaluated at p = Pres- 

Equation |T^ is a non-linear equation in The righthand side of this equation 
depends on through p^es directly, through the dependence of X on p^es and through 
the number densities. In order to solve this equation, we need to write the chemical 
potentials in terms of Ly^. Now, for each temperature T, the neutrino asymmetry is 
created at the neutrino momentum Pres- However, for temperatures greater than about 
1 MeI/[^ the effect of the weak interactions is to quickly thermalise the neutrino 
momentum distributions. This means that the neutrino asymmetry is approximately 
distributed throughout the neutrino momentum spectrum via chemical potentials for 
the neutrinos and anti-neutrinos. In general. 


= 


x^dx 


x’^dx 


4^ ^3^ Jq gx+Mz/ 4^ ^3^ Jq 


( 20 ) 


where hi= Pi/T and i = Expanding Eq. (pO]) , we hnd that 


Ly ~ 


24C(3) L 


- hp) - 6(/i^ - Py) ln2 + {p^ - p^) 


( 21 ) 


which is an exact equation for Py= — Pp, otherwise it holds to a good approximation 
provided that Py^p ~ 1. For T ~ (where ~ 3 MeH and ~ 5 MeH are 
the chemical decoupling temperatures) py^ ~ ~Pua because processes such as Ua + 

r\j 

Pa + e“ are rapid enough to make Py + Pp Pf.+ + he- — 0. However, for 

1 MeV ~ T ~ Tdec^ weak interactions are rapid enough to approximately thermalise 
the neutrino momentum distributions, but not rapid enough to keep the neutrinos in 
chemical equilibrium. In this case, the value of Py is approximately frozen at T ~ 
while the anti-neutrino chemical potential Pp continues increasing until T ~ 1 MeV. 


We also need to specify the initial condition in order to solve Eq. dlSf) . To do so 
we need to know the value of Ly^ at some temperature Tj < Tc at which MSW tran¬ 
sitions are already dominant. This Ly^ value can be obtained by solving the exact 
quantum kinetic equations, based on the density matrix, which incorporate both col¬ 
lision and oscillation effects. Fortunately, it turns out that the subsequent evolution 
of lepton number is reasonably insensitive to what temperature Tj is chosen as the 
initial temperature for Eq.(p!8[), provided that T, is chosen during the epoch after Tc 


for which Ly^ 1[0- When the asymmetry Ly^ -C 1, Eq.([l^) can be simplihed to 


dLyjdT ~ -ALy 


JT which means that Ly^T^ 


is approximately constant. This means 


that Pres/T is also approximately constant given that Ly^ is related to the resonance 
momentum p^es by Eq. ([1T|) . As we will discuss in section V, a numerical solution of the 
quantum kinetic equations shows that Pres/T is generally in the range 


0.2 


~ Pr 


,IT ~ 0.8, 


( 22 ) 


for T values around Tj = Tc/2 when the oscillation parameters have been chosen to 
lie in the parameter space of interest, which turns out to be \5m\g\ ~ 10“^ eV"^ and 
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sin^26'o ~ 5 X /\5rn^^g\Y/^. (This lower bound for the mixing angle ensures 

that a suitably large asymmetry is created at Tc 0. The result of Eq.(^) can also be 
gleaned from the static approximation based results of Ref. |^.) We will from now on 
use a value of about Tc/2 for Tj. 

Before presenting the results of a numerical solution of Eq. (|T^ for the hnal asymme¬ 
try , it is interesting to note that an approximately correct answer is easily obtained 
from the following argument. As T falls below Tc/2, the asymmetry keeps increasing. 
This eventually forces the rate of change of to decrease substantially. Recall that 
dLy^/dT is proportional to the how quickly the resonance momentum pres moves as per 
Eq.(p^. When is large, Pres/T must move to large values in order to create lepton 
number. Eventually, Pres/T —> oo, and all of the anti-neutrinos which have passed 
through the resonance have been converted into sterile neutrinos. Thus, assuming that 
the initial number of sterile neutrinos is negligible and also neglecting the modihcation 
of the distribution due to the chemical potential, we expect that the hnal value of the 
lepton number, , to be given roughly by 


Ll^ 1 x^dx 3 

h 4C(3) 4i„/T 1-h 8’ 


(23) 


where h = (note that h ~ 1 for T ~ rUg — 0.5 MeV) and Pin/T is the value 

of Pres/T [and is in the range of Eq.(|^)] at T ~ Tc/2. It is interesting that the hnal 
asymmetry is approximately independent of Pin/T and hence also of This is 

because Pin/T from Eq. (|^) is always small. 

Actually, the hnal value of the lepton number is somewhat less than 3/8 = 0.375 if 
it is created when T ~ 1 MeV. This is because the number density of anti-neutrinos 
is continually reduced as the lepton number is thermally distributed via the chemical 
potential. Thus, depends on the temperature at which Ly^ becomes large (10“^ 
roughly), and thus on \drn/^^\. Numerically solving Eq. (|I8|) , assuming that the initial 
number of sterile neutrinos is negligible, we found that the hnal value of the lepton 
number is[P^ 

Li^/h ^ 0.35 for \5mlf/eV‘^ ~ 3, 

L/^/h ~ 0.23 for 3 ~ ~ 3000, 

Lljh ~ 0.29 for \6mlf/eV^ ~ 3000. (24) 

In numerically solving Eq. (p!8D we start the evolution at T ~ Tc/2 with Pin/T in the 
range of Eq.(|2^) and with a corresponding Ly^ obtained through Eq. (PR|) . We found 
that is approximately independent of the initial value of Ly^ for Pin/T in this range. 

The temperature where the hnal neutrino asymmetry is reached is approximately, 

T/~ 0.5 (|W|/eV^)^^^ MeV. (25) 

This result can be obtained analytically by using the resonance relation Eq.(^) with 
Ly^ ~ and Pres/T ~ 6 (since is not reached until Pres/T 3> 1 and we take 
Pres/T rsj 6 for dehniteness) . 
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Equation is an approximation based on the neglect of collisions and the as¬ 
sumption of complete MSW conversion. By numerically integrating the exact quantum 
kinetic equations]^, we have checked that Eq.(P!B[) does indeed accurately describe the 
evolution of the neutrino asymmetry in the range 1 MeV ~ T ~ Tc/2. We will 
discuss this and provide an illustrative example in section V. 

As preparation for the application of the above formalism to BBN, we need to 
discuss how an asymmetry in Vf, can be generated in the context of an overall neutrino 
mixing scenario. There are two generic ways of producing a nonzero Ly^. First, — Vg 
oscillations can generate Ly^ directly. Alternatively, — (and/or v^ — Vs) oscillations 
can generate a large Ly^ (and/or Ly^) some of which is then transfered to Ly^ by Pr — 
(and/or — z/g) oscillations. 

The direct way of generating Ly^ is only possible in special circumstances. Either 
l^rUg^l S> \6m‘^g\, or Ug has signihcant mixing with z/g only. Only in these cir¬ 

cumstances can we consider z/g — Ug oscillations in isolation. For this case we have 
estimated the effects of the neutrino asymmetry on BBN by writing a nucleosynthesis 
code. We hnd that —1.8 ~ 5Nyf^ ~ —0.1 requires a \5w?^g\ in the range 0.5 — 7 eE^. For 
I^TTig^l ~ 0.5 lepton number is created too late to signihcantly affect BBN, while 
for l^rUg^l ~ 7 lepton number is created so early that it leads to 5Ny ~ —1.8 and 
thus appears to be too great a modihcation of BBN to be consistent with the observa¬ 
tions. Note however that for sin^ 26q large enough, the sterile neutrino can be excited at 
temperatures before signihcant lepton number is generated (which for \SmU ~ 1 eV is 
T ~ 13 MeV). This can lead to an increase in the energy density which can (partially) 
compensate for a large positive electron lepton number. 

While the above direct way of generating Ly^ is a possibility, we believe that a 
more interesting possibility is that Ly^ is generated indirectly. As we will show, this 
mechanism gives 6Ny ~ —0.5 (assuming Ly^ > 0) for a wide range of parameters. This 
mechanism is also the only possibility if \Sm^g.\ S> ^ assuming 

that Pg mixes with all three ordinary neutrinos. 


IV. An example with four neutrinos 

Consider the system comprising Pr,p^,Pe,i^s- An experimental motivation for the 
sterile neutrino comes from the current neutrino anomalies. There are several ways 
in which the sterile neutrino can help solve these problems. For example, the solar 
neutrino problem can be solved if Sml^/eV'^ ~ 10“® and sin^ 20o — 10“^ (small angle 
MSW solution)P or if 10“^ A l^rUg^l/eV^ ~ 10“^° andsin^26*o — 1 (maximal oscillation 
solution)!^. Alternatively, p^ — Pg oscillations can solve the atmospheric neutrino 
problem if ~ 10“^ eV^ and sin^ 26*o — 1 

We will assume that niy^ ^ niy ^, rriy ^, niy ^, which means that 

\5M‘^\ = |5m/g| ~ ~ > \5mlX |(5mjg|. (26) 

With the above assumption, p^ — Pg oscillations initially create signihcant Ly^ at the 
temperature T = T^. — 16(|5M^|/eV^)^/®. As before, we will assume that the sign of 


lA 
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is positive]^. The effect of Vr — v^, and oscillations is to generate and 

Ly^ in such a way that = Ly^ — Ly^ —*> 0 and = Ly^ — Ly^ —*> 0, 

respectively [^, |^. [Note that if Ly^ > 0, then the MSW resonances for Vj. — v^. and 

Ur — ^fi oscillations occur for anti-neutrinos (given also our assumption that niy^ > rriy^^^) 
and so the signs of Ly^ and Ly^ are also positive]. However, the rate of change of lepton 
number due to collisions, the dominant process at higher T, is typically too small to 
efficiently generate Ly^ from L.Jli. But, as Ly^ becomes large at lower T, lepton 
number can be efficiently transferred by MSW transitions. (When Ly^ 1, MSW 
transitions cannot efficiently create Ly^ because Ny^ — Ny^ ~ 0 and MSW transitions 
only interchange Vr with z/g without changing their overall number density). The rate 
of change of lepton number due to Va — oscillations is simply given by the difference 
in rates for which z/q, anti-neutrinos and anti-neutrinos pass through the resonance 
(assuming that sin^ 26q 1). We need to consider the three resonances, z/,- — 9^- — Pg 

and P,- — P^, for our system. We denote the resonance momenta of these resonances 
by Pi, P 2 and pa, respectively. The rate of change of the lepton numbers due to MSW 
transitions is governed approximately by the differential equations. 


dLy^ .^1 d{pi/T) I -y I d{p2/T) I y \ d{p3 /T) i 

PUP — ^l| PUP \ -^21 pup: \ Aal I, 

W I d{p 3 /T) I dLy^ = J-X I I 

dT “Tv^ai pT I’ dT prp I, 


(27) 


where 







ZVo.). .’f2 



. .’fa 



(28) 


and the Xj are evaluated at p = pj {i = 1,2,3). Note that Xj depends on T through 
the ratio Pi/T and through the dependence of the various chemical potentials on T. 
Observe that 


d{PrlT) ^ d{p,/T) d{pi/T) dLy^ d{pi/T) dLy^ djp^/T) dLy^ 

dT dT dLy^ dT ^ dLy^ dT dLy^ dT ’ 


with 


diPi/T) ^ n d(pi/T) ^ n d(pi/T) ^ -2(pi/T) 
dLy^ ^ ^ dLy^ W) ’ 

d(P2/T) d(p2/T) -(P2/T) d(p2/T) p 

dLy^ dLy^ L(^)-L(0 ’ dLy^ 

d(P3/T) ^ d(p3/T) ^ -(p 3/T) d(p3/T) ^ r. 

dLy^ dLyy^ dLy^ 

ajpi/T) ^ -ipi 
dT ~T^- 

By the symmetry of the problem, Ly^ = Ly^, p 2 = pa and dLy^/dT 
Using this simplihcation, we hnd that 

dLy^ dLy^ y4. dLy^ OC [3 dLy^ 


(30) 

dUJdTB- 

(31) 
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where A = jyi + aS, B = yiy2 — [ 36 , with 


PI 


P2 

R ox. 

Pi 


P2 

T2 

oyv2 

T'^ 

1 H 

tURR 


T{LR)-L(R)_ 


7 — 4X2 


P2_ 

T2 


2/1 — 1 + 2X1 


E3 


TLt^) 


+ 2 X 2 


, 5 = +X 2 

P2 

.r(L(^)-L('=)). 


P2 

, 1/2 = 1 + X 


P2 

2 [r(L(^)-L('^)) 


( 32 ) 


In deriving this equation we have assumed that d{pi/T)/dT < 0 and d{p2/T)/dT < 0. 
Observe that X2d{p2/T)/dT = —A/B and thus for selfconsistency Eq.(PI|) is only valid 
provided that A/B <Q (given that X2 < 0). If d{p2/T)/dT > 0 then Eq. (^ becomes 


dLy^ dLi/^ A dL^^ 01 jd dL^^ 
~dT^~d/F^^ ~d/F ^ J^~dT' 


(33) 


~ ~ /NJ 

where A, B, a, /?, yi have the same form as A, B, a, /?, yi except that X 2 —*• —X 2 . In this 
case, X 2 d{p 2 /T)/dT =A / B- It follows that Eq. (pSD is only selfconsistent provided that 
A / B< 0. Observe that d{p 2 /T)/dT must be continuous which means that d{p 2 /T)/dT 
only changes sign when A changes sign and Eq. (|3T|) maps onto Eq.(^^ continuously 
because A= —A (and thus A= A at the point where A = 0). If d{pi/T)/dT changes 
sign at some point pi/T = q then we must make the replacement Xi ~ 0 for pi/T < q 
(assuming that initially d{pi/T)/dT < 0) since the previous MSW transitions have 
populated Ug for pi/T < q. 

In solving Eq. (|3TD , we will assume that the initial number of sterile neutrinos can be 
neglected (this will be valid for a wide range of parameters as will be discussed later). 
We start the evolution of Eq. (|3T|) when T ~ Tc/2 (with T^. given by Eq.(|^ for Vj. — Vg 
oscillations). There is a range of values of at this point which is related to the range 
oipres/T [Eq.(^)] through Eq. (pTf) . Performing the numerical integration, we hud that 
the hnal electron neutrino asymmetry is [^] 

Lljh ~ 2.0 X 10-2 for 10 ~ \5M^\/eV^ ~ 3000 

Lljh ~ 1.7 X 10-2 for |5M2|/eI/2 ~ 3000 (34) 

and recall that h = Tf^/T^. We found that is approximately independent of Pres/T 
for Pres/T in the range given by Eq. (p2D . We also found numerically that Ll^ is ap¬ 
proximately independent of the initial value of Ly^ (at T ~ 77/2), so long as Ly^ ~ Ly^ 
at this temperature (which should be valid since efficient generation of Ly^ does not 
occur until much lower temperatures where Ly^ has become very large). In addition, 
we found that is independent of the precise value of the initial temperature (so long 
as the initial temperature is less than T^ and Ly^ -C 1 at this temperature). The reason 
for this independence is simply due to the fact that signihcant generation of Ly^ cannot 
occur until Ly^ becomes large (~ 10-2). The hnal asymmetry is also independent of 
sin2 29q so long as sin2 29q -C 1 for aforementioned reasons. Finally, and perhaps of most 
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interest, we find that Ll^ is almost independent of |hM^| so long as |5M^| rsj 3 eV^. For 
|(5M^| ~ 3 eV^, does not become large until T ~ 1 MeV. For temperatures in this 
range, the effect of cannot be described in terms of chemical potentials because the 
weak interactions are too weak to thermalise the neutrino distribution. For this reason, 
Ll^ should be much smaller since the Ur — resonance (which occurs at a momentum 
which is always greater than the — Vg resonance) simply interchanges almost equal 
numbers of i^t-’s and i^e’s. 

We now apply the above analysis to BBN. Recall that the neutrino oscillations affect 
in two ways. First, the creation of and the related modihcation of the neutrino 
momentum distributions directly affects the nuclear reaction rates which determine the 
neutron/proton ratio. Second, the oscillations can modify the energy density of the 
Universe by the excitation of the sterile neutrino and the modihcation of the neutrino 
momentum distributions due to chemical potentials. We hrst discuss the energy density 
question. 

For T > Tc the Vr — i^s oscillations can excite the sterile neutrino (and anti-neutrino). 
In Ref.p, a detailed study was done which found that pus!Pv ~ 0-6 provided that 


sin^ 200 ~ 4 X 10"® 


r eU2 


1/2 


|(5M2 


(35) 


Furthermore, PvJPv 
Ps/Pu ~ 0.1 for 


0 very quickly as sin^ 20o —>-0. In particular, we found that 


sin^ 200 ~ 5 X 10"® 


r eU2 


1/2 


|hM2 


(36) 


Note that after the lepton number is created, the oscillations no longer excite signihcant 
numbers of sterile neutrinos until ~ 10“^. At this point the Vj- — Pg oscillations 
(recall that we are assuming that Ly^ > 0) transfer Vt —>■ Pg. The effect of these 
oscillations on the overall energy density depends on the temperature where Ly^ ~ 10“^ 
occurs, which in turn depends on There are essentially three regions to consider, 

10 ~ \5M^\/eV^ ~ 3000, \5M‘^\/eV‘^ ~ 10 and \5M‘^\/eV‘^ ~ 3000. 

For 10 ~ |hM2|/eU2 3000, we have numerically calculated the hnal number 

and mean energies of Pr,Ps,Pe,i^ii (the number and energy densities of the neutrinos 
are approximately unchanged in this region). Normalizing the number density to the 
number of neutrinos when py = 0, no = |C(3)T^/7r2, 


—^ = —- ~ 0.95, —^ ~ 0.44, —^ ~ 0.66. (37) 

no no no no 


Note that the total number is approximately unchanged (i.e. 0.95 x 2-1-0.44-1-0.66 ~ 3). 
We hnd the hnal mean energy for the Pg, {Eg), to be slightly less than than the mean 
energy for a Fermi-Dirac distribution with py = 0, (Ep^) ~ 3.15T [{Eg)/{Epp) ~ 0.88]. 
For this reason there is a small overall change in energy density, equivalent to about 
6N/E ~ —0.05. For |hM2| ~ 3000 eU^, Ll^ is reached for T ~ and so py^ ~ 
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—jj^u-r- III this case, there is an additional contribution to the energy density coming 
from the Uj. neutrinos due to the negative chemical potential In this case we hnd 
that the overall change in the energy density is considerably larger and equivalent to 
~ 0.4. Finally, for |5M^| ~ 10 eld^, the change in the energy density quickly 
becomes completely negligible because the weak interactions are unable to thermalise 
the neutrino distributions. The oscillations simply transfer to Vg and the total number 
and energy density remain unchanged. 

We now turn to the effect of and the corresponding modihcation of the momen¬ 
tum distributions on through nuclear reaction rates. For 10 ~ \5M‘^\/eV'^ ~ 3000, 
the distribution of Ll^ can be approximately described by chemical potentials hp ~ 0.06 

and hi/ — 0. For ~ 3000 eld^, the lepton number is created above the chemical 

decoupling temperature. In this case the distribution of Ll^ can be approximately 

described by chemical potentials hp — 0.025 and hi/ — —0.025. We hnd that for 
\5M‘^\ ~ 10 eld^, L = Ll^ is reached for T ~ 1.5 MeV. Thus to a good approxi¬ 
mation the chemical potentials hi/e,hpe are approximately constant during the nucle¬ 
osynthesis era. Using our BBN code, we hnd that the modihcation of Yp due to the 
chemical potentials is 6Yp ~ —0.005 for hp — 0.06, hi/ — 0 and 6Yp ~ —0.006 for 
hp — 0.025, hi/ — —0.025. From Eq.(|^), this translates into a reduction of the ehective 
number of neutrino degrees of freedom during nucleosynthesis. Including the ehects of 
the change in energy density discussed earlier, we hnd that 

^ _o.5 for 10 ~ \5M^\/eV^ ~ 3000 
SN^ff ~ _o.l for \5M^\/eV^ ~ 3000. (38) 


For this result, we have considered the case of negligible excitation of sterile neutrinos 
for temperatures above Tc, that is, Eq.(|3BD has been assumed. Note that if we had 
assumed that was negative instead of positive, then the sign of is also negative 
and the change in Yp due to the asymmetry is opposite in sign as well. This leads to 
5 N^ff ~ +0.4 (0.9) for 10 ~ \5M^\/eV^ ~ 3000 {\5M‘^\/eV^ ~ 3000). 

In our analysis we have neglected the ehects of the — Pg) ~ ~ 

Vg oscillations. It is usually possible to neglect these oscillations if \5m?\ -C \5M‘^\ 
because the lepton number created by Pr — oscillations is large enough to suppress 
the oscillations which have much smaller Srn^. Of course, in some circumstances these 
oscillations cannot be neglected. For example, in Ref. [^, we showed that the ehects of 
maximal — Ug oscillations with ~ 10“^ eU^ (as suggested by the atmospheric 

neutrino anomaly ||24||) can only be neglected if ~ 30 eU^. Interestingly, this 

parameter space overlaps considerably with the parameter space where ~ —0.5, 

according to Eq.(^). Note that this parameter space is also suggested if the tan 
neutrino is a signihcant component of dark matter. 
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V. Evolution of lepton number from the exact quantum kinetic eqnations 


In this section we study the evolution of the neutrino asymmetry by numerically 
integrating the exact quantum kinetic equations!^. This formalism allows a near 
exact calculation to be performed which is valid at both high and low temperatures. 
As we have discussed, for high temperatures T ~ Tc the evolution of lepton number 
is dominated by collisions (assuming \6rn?‘\ ~ 10“^ eV"^) while at lower temperatures, 
the evolution of lepton number is dominated by oscillations between collisions (MSW 
effect). 

The system of an active neutrino oscillating with a sterile neutrino can be described 
by a density matrix |T8|, 28[| . Below we very briefly outline this formalism. The den¬ 
sity matrices describing an ordinary neutrino of momentum p oscillating with a sterile 
neutrino are given by 


pM = \Po(p)\l + P(p) ■ ff], Pp(p) = i^Po(p)[l + P(p) ■ ff], 


(39) 


where P(p) = Px(p)£ + P^(p)p + Pz(p)z. (It will be understood throughout this section 
that the density matrices and the quantities Pi also depend on time t or, equivalently, 
temperature T.) The number distributions of z/q, and Ug ci-re given by 


N.. = i/>o(p)|l + PMKl N,. = ifo(p)|i - P.(p)]Kl 


where 


= 


p 


27r2 1 + 


exp 


( 

V T 


y 


(40) 


(41) 


is the equilibrium number distribution. Note that there are analogous equations for the 


anti-neutrinos (with P(p) - 
governed by the equations 


P(p) and Pq —> Pq). The evolution of Po(p) and P(p) are 


Wt 


d 


P(p) = V(p) X P(p) -h [1 - Pz{p)] 




InPo(p) 


d 


^p«(p) - 


[D{p) + ^ lnPo(p)][Px(p)x + Py{p)y], 
Rip). 


The quantity V [p) is given by 

V(p) = /3(p)x + A(p)z, 
where Pip) and A(p) are defined by 

6rrP 


Sm ^ 

f^ip) = ^ sin 200, A(p) = 

Ip 


2p 


cos 200 — Kp) a 


(42) 


(43) 


(44) 
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in which the +(—) sign corresponds to neutrino (anti-neutrino) oscillations. The di¬ 
mensionless variables a{p) and b{p) contain the matter effects and are given in Eq.(|^). 
The quantity D{p) is the quantum damping parameter resulting from the collisions of 
the neutrino with the background. According to Ref.[^, the damping parameter is 
half of the total collision frequency, i.e. D{p) = r^^{p)/2. Finally, note that in Eq.(^^ 
the function R{p) is related to T^^{p) and its specific definition is given in Ref.[18|. For 


temperatures above 1 MeV, we can make the useful approximation of setting 

and Np^ = N^. This means that P^ip) = 2/[l + Pzip)], Poip) = 2/[l + Pzip)] and 

consequently 


dPojp) _ -2 ___ 

dt [l + Pzip)]"^ dt ’ dt [1 + Pz{p)]'^ dt 


dP^{p) dPo{p) 


dPz{p) 


(45) 


For the numerical work, the continuous variable p/T is replaced by a finite set of 
momenta Xn = Pn/T (with n = The variables P(p) and Po{p) are replaced 

by the set of N variables P(a:n) and Po{xn)- The evolution of each of these variables 
is governed by Eqs.(|4^), where for each value of n, the variables V(p) and D{p) are 
replaced by V (xn) and D{xn)- Thus, the oscillations of the neutrinos and anti-neutrinos 
can be described by 8N simultaneous differential equations. 

The rate of change of lepton number is given by 



1 

{Uu^ -npj 

1 

1 

{n^s - np,) 

dt 

1 

a- 

fTH- 1 

_1 

1 

a- 

fTH- 1 

-1 

e 

_1 


Thus, using Eq. (|40|) , 


cLLu^ \ d 1 

dt 2 dt 


Foil-P z)N^pI-F oil 


P.)NZ dp 


Taking the time differentiation inside the integral we find that 


^ ^ 1 /■ / a[-Po(l - ft)] NZ _ c)|Po(l-P.)] iVa\ 

dt 2 J \ dt Tiry dt rij ) 


(46) 


( 47 ) 


(48) 


where we have used the result that N^^dp/n^ is approximately independent of t. Ex¬ 
panding this equation using Eq. (^5[) , we find 


dLij 


dt 



2 ap, 

- [1 + PX St '11 + p,f dt , 


dp 


(49) 


30 


Equations ^ and can be numerically integrated to obtain the evolution of L, 

. We illustrate this with an example. For definitness we will consider the Ur, Vs 
system. In Figure 1, we take = —10 and sin^ 200 = 10“® (we set r/ = 4 x 10“^° and 
took = 0 initially]^). The result of numerically integrating Eq. (^2D and Eq.(^) is 
shown in the figure by the dashed-dotted line. Also shown in Figure 1 (dashed line) is the 
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“static approximation” (Eqs.(94) and (93) of Ref. |^). As discussed in Ref. [Q, the static 
approximation assumes that the system is sufficiently smooth and that the dominant 
contribution to the rate of change of lepton number is collisions. As shown in Figure 1, 
the static approximation is a good approximation at high temperatures. However, as 
discussed in Ref. p, the static approximation does not include the MSW effect which is 
the dominant physical process at low temperatures. As expected the MSW effect keeps 
growing like Ly^ ~ 1/T^ for much lower temperatures. We have also checked our 
simplihed Equation (p!8D for the evolution of lepton number due to MSW transitions. 
We started the evolution of this equation at T = Tc/2 ~ 13.5 MeV with the value of 
Ly^ at this point obtained from the quantum kinetic equations of Ly^ ~ 2.92 x 10“®. 
The subsequent evolution of Ly^ obtained from numerically integrating Eq.(^) is given 
in Figure 1 by the solid line. As the hgure shows, Eq.(]T^) is a very good approximation 
for the evolution of the neutrino asymmetry at low temperatures. This provides a useful 
check of the validity of the approximate approach used in section IV for the z/g, z/^, Vg 
four flavour system. 

It is instructive to examine the evolution of the neutrino and anti-neutrino resonance 
momenta. Recall that the resonance for neutrinos occurs when b — a = cos 200 while 
the resonance for anti-neutrinos occurs when h + a = cos20o- Let us write 

b = Aip^, a = X 2 P, (50) 


where Ai and A 2 are independent of p and can be obtained from Eq.(^. Note that 
Ai, A 2 > 0 given that Sm^g < 0 and assuming > 0. Solving the resonance conditions 
b ±a = cos 200 , we hnd that the resonance momenta satisfy 


^2~^\/ AqH-4Ai cos 20q p 

Pres =- 2 Ai - neutrinos^ 

—A2H-4Ai cos 20q p . . 

Pres =- ^ 2 Ai - “ neutrinos. 


(51) 


In Figure 2, we have plotted the evolution of the resonance momenta for the neutrinos 
and anti-neutrinos. As this example illustrates, the neutrino resonance momentum 
moves to very high values as T ~ Tg while the anti-neutrino resonance momentum 
moves to very low values (which in this example is Pres/T ~ 0.6 for T ~ Tc/2). We 
have found that this behaviour is quite general, with the anti-neutrino resonance Pres/T 
in the range Eq. (|^) at T = Tc/2 as sin^ 20o and 5m^ are varied. 


VI Conclusion 


In summary, we have extended previous work on neutrino oscillation generated 
lepton number in the early Universe by studying the evolution of lepton number at low 
temperatures where the MSW effect is important. We applied this work to examine 
the implications of the neutrino asymmetry for BBN in two illustrative models. In the 
hrst model, electron neutrino asymmetry was created directly by z/g — Vs oscillations 
while in the second model the electron neutrino asymmetry was created indirectly by 
the reprocessing of a tau neutrino asymmetry. 
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One result of this study is that the naive conclusion that sterile neutrinos only in¬ 
crease the effective number of neutrino species during the nucleosynthesis era is 

actually wrong. Neutrino asymmetries generated by neutrino oscillations can naturally 
lead to a decrease in . Furthermore in the case where the electron neutrino asym¬ 
metry is transferred from the tau or mu neutrino asymmetries, the electron neutrino 
asymmetry is approximately independent of and sin^ 26^0 for a wide range of pa¬ 
rameters. This leads to a prediction of ~ —0.5 if the asymmetry is positive for 

an interesting class of models. Remarkably, this prediction is supported by some recent 
observations which actually suggest < 3. 
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Figure Captions 


Figure 1. The evolution of the — Ug oscillation generated lepton number asymmetry, 
We have taken by way of example, the parameter choice 6m‘^ = —10 and 
sin^ 200 = 10“®. The dashed-dotted line is the result of the numerical integration of 
the quantum kinetic equations [Eq.(^) and Eq.(^^]. The solid line is the result from 
the numerical integration of Eq.(|l^) while the dashed line is the static approximation 
developed in Ref. [Q. 

Figure 2. The evolution of the neutrino (dashed line) and anti-neutrino (solid line) 
resonance momenta for the example of Figure 1. 
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